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We discuss the possibility of finding n-type nematic order in the vicinity of a FM phase 
in the multiple-spin exchange model on the triangular lattice. We study this problem both 
from 5 = 00 (classical) and 5 = 1/2 (quantum) limits. 

The possibility of realizing a quantum spin liquid in a frustrated magnet with 
antiferromagnetic interactions has been widely discussed. Here we consider the na- 
ture of a quantum spin liquid appearing in the vicinity of a ferromagnetic phase. 
Our motivation is the gapless spin liquid observed in two-dimensional (2D) solid ^He 
on graphite/^' ^-^ very close to a ferromagnetic phase. This system is believed to be 
well described by a multiple-spin exchange (MSE) model on the triangular lattice: 

n = jY,P2 + K ^(P4 + p^') -hY, si, (1) 

N.N. p i 

where Pn denotes the cyclic permutation operator of n spins. The first summation 
runs over all nearest neighbor pairs and the second one all four-site minimal diamond 
plaquettes. In 2D solid '^He, the effective two-spin exchange coupling J < is 
ferromagnetic (FM) while the four-spin cyclic exchange K > is antiferromagnetic 
(AF). An exact-diagonalization study of finite-size systems^^ concluded that the FM 
phase for small K evolves into a state with zero total spin for large values of K and 
that, if K is quite large, the ground state is spin liquid with no long-range magnetic 
order and a spin gap of order \J\. How such a spin liquid might form out of a FM 
state, and what the nature of the spin liquid should be in the vicinity of the FM 
phase, remain open problems. 

Here, we present evidence that the first instability of the FM state with increas- 
ing K is against a gapless spin liquid state with nematic order. We consider the 
n-type spin nematic order parameter, defined by 

0^^{n,r,) = i(5r5f + S^Sf) - ^6»P{S, • S,). (2) 

This type of spin nematic order has been seen in the S = 1 bilinear-biquadratic 
model. However, to the best of our knowledge, it has never before been found in a 
spin S = 1/2 frustrated system. 

Let us consider first the classical limit of the MSE model on the triangular 
lattice. A mean-field argument and its extension^) showed that a highly degenerate 
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Fig. 1. (a) Temperature dependence of specific heat in the classical MSB model for J = —1, K = 
1/3, H = with peak corresponding to the onset of nematic correlations. Results are shown 
for systems of A'^ = 432, 768 and 1200 spins, (b) Decay of spin (circle) and nematic (triangle) 
correlations at T = 0.05 for the classical MSE model with J — —1, K = 1/3, H = 0, and 
A'^ — 1200. All temperatures are measured in units of \ J\. 



phase appears in the wide parameter range — 1 < K/J < —1/4. This set includes 
all parameters relevant to 2D solid ^He in the low density limit, and extends up to 
the classical boundary of the FM phase at K/J = —1/4. Among the many possible 
classical ground states, collinear states are most likely to be favored by quantum or 
thermal fluctuations and the two most characteristic collinear ground states are the 
UUUD and UDDD states with exactly half the saturation magnetization. However, 
because it costs no energy to add a straight domain wall between these states, at 
any finite temperature the system will gain entropy by breaking up into stripe-like 
UUUD and UDDD domains. For this reason we do not expect any long range spin- 
spin correlations to form in this system, even in the limit T — > 0. 

However we do expect that the spins will remain collinear, and that this will 
lead to long ranged spin nematic correlations of the form 0"^(rj,rj), above. To 
test this argument, we performed MC simulations of the classical model and found 
that nematic correlations are enhanced at low temperatures, while spin correlations 
decay quite rapidly. (See Fig. [TJ) The onset of nematic correlations is signaled by a 
peak in the specific heat whose height remains finite as the system size is increased. 
Of course, true long ranged nematic order cannot be achieved in 2D at any finite 
temperature because the director (here, the axis of collinearity) breaks a continuous 
symmetry. 

We now turn to the experimentally relevant quantum limit, i.e. the 5=1/2 MSE 
model, and consider the nature of the first instability of the FM as K is increased. 
In the absence of magnetic field, the entire one magnon spectrum 



e{k) = -2( J + AK){3 - cos(fc • ei) - cos(fc • eg) - cos[fc • (ei - 62)]} (3) 
vanishes at the classical phase boundary K/J = —1/4, reflecting the extensive de- 



Nematic order in the multiple- spin exchange model 



3 




Fig. 2. Saturation fields: the parameter K/\J\ vs. the one- and two-magnon saturation fields He 
and Hc2- In the shaded region, the bi-magnon bound states have a negative energy relative to 
the saturated FM state, while individual magnons have a positive energy with a finite positive 
stiffness. Thus the partially polarized spin state in the shade should show BEC of bound bi- 
magnons. 

generacy of the classical disordered phase beyond it. For large values of applied 
magnetic field, this degeneracy is lifted and the first instability of the saturated FM 
state is against a three-sublattice canted AF state. However for smaller applied 
fields, the system can gain kinetic energy by allowing pairs of fiipped spins to move 
together through repeated use of the cyclic exchange process. This prompts us to 
consider the following trial bi-magnon bound state: 

a{k)\l{k)) + b{k)\2{k)) (4) 

with 

\m) 
\m) 

Despite the fact that the static interaction between pairs of magnons is repulsive, 
for K/J < —4/17, there exist trial solutions with a negative energy relative to the 
saturated FM state, i.e. bi~magnon bound states below the FM continuum. Thus 
the first instability of the saturated FM state in zero field is against the spontaneous 
creation of bound bi-magnon pairs. Extending this calculation to finite magnetic 
field, we obtain the phase diagram shown in Fig. [2j 

In order to better understand the unsaturated phase bordering on the FM phase 
let us consider the dynamics of bound bi-magnons. We know that pairs of fiipped 
spins have repulsive interactions and a definite (positive) hopping amplitude. From 
this we can write down an effective Hamiltonian for bound bi-magnons, which has 
a form similar to an AF XXZ model. We therefore anticipate that the density of 
bound bi-magnons increases continuously upon decreasing the magnetic field from 
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the critical value Hc2, and that there is a second order quantum phase transition 
from the saturated to the unsaturated state. 

To be certain that bound bi-magnons arc the elementary particles of the un- 
saturated state, we need to rule out the existence of larger magnon bound states. 
Cooperative motion of more than three flipped spins can occur in high order pertur- 
bation theory in K. However the repulsive interaction between magnons is first order 
in K, so we do not, in general, expect larger bound states to be stable. From these 
arguments we conclude that the unsaturated states close to the saturation field are 
characterized by the Bose-Einstein condensation (BEC) of a low density of bound 
bi-magnons. 

Finally we consider the nature of this bound bi-magnon condensate. The magnon 
pairing operator is given by 

Op^ir{ri,rj) = s-sJ. (5) 
This is related to the nematic operators through 

Opair = srsj = ow - - 2io%y, (6) 

where we have omitted site indices. The real part of the pairing operators corre- 
sponds to the nematic operator Off — Oj^ and the imaginary part to —20^. These 
two components give nematic order parameters for the systems under the magnetic 
field. Thus, the BEC of bound bi-magnons is equivalent to the emergence of long- 
ranged spin nematic order. 

Our conclusion is that the MSE model on the triangular lattice in the vicinity of 
the saturated FM phase shows long ranged spin nematic order in one of two channels 

, hm ([O^^(0,ei) - O^f(0,ei)][O--(n,r, + ei) - 02'^(r,,n + ei)]) 

ci{l - m/msat) (7) 

or 

hm 4(0^2^(0,61)0^^(^,^ + 61)) - ci{l-m/msat), (8) 

>oo 

where m (nisat) denotes the (saturated) magnetization and ci a finite constant. In 
these states, spin correlations decay exponentially even at zero temperature, making 
them as true spin-liquids. The rate at which the spin correlations decay is set by 
the binding energy of bound bi-magnon pairs, which vanishes at the boundary with 
the saturated FM phase. 

Correlated hopping processes^) lead to a similar BEC of bound bi-magnons in 
the low-magnetization regime of the 2D Shastry-Sutherland model^-* . This can also 
be regarded as a spin nematic state. However, nematic order may be difficult to 
observe in the "Shastry-Sutherland" compound SrCu2(B03)2 because of the pres- 
ence of Dzyaloshinsky-Moriya interactions which break the bi-magnon pairs. The 
spin interactions in solid ^He, on the other hand, are purely isotropic in nature. 
We therefore expect that nematic order is experimentally accessible for a range of 
densities of 2D solid ^He in applied magnetic field. 
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We will return to this issue, and present further analytic and numerical results 
for nematic order neighboring ferromagnetism in the MSE model, elsewhere.^^ 
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